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DISCOVERY OF A WANDERING RADIO JET BASE AFTER A LARGE X-RAY FLARE IN THE BLAZAR 

MARKARIAN 421 

K. Niinuma\ M. Kino^, A. Doi^, K. Hada'^, H. Nagai® and S. Koyama® 

Not to appear in Nonlearned J., ^5. 

ABSTRAGT 

We investigate the location of the radio jet bases (“radio cores”) of blazars in radio images, and 
their stationarity by means of dense very long baseline interferometry (VLBI) observations. In order 
to measure the position of a radio core, we conducted 12 epoch astrometric observation of the blazar 
Markarian 421 with the VLBI Exploration of Radio Astrometry at 22 GHz immediately after a large 
X-ray flare, which occurred in the middle of 2011 September. For the first time, we find that the radio 
core is not stationary but rather changes its location toward 0.5 mas downstream. This angular scale 
corresponds to the de-projected length of a scale of 10^ Schwarzschild radii (J?s) at the distance of 
Markarian 421. This radio-core wandering may be a new type of manifestation associated with the 
phenomena of large X-ray flares. 

Subject headings: galaxies: active — galaxies: jets — BL Lacertae objects: individual (Markarian 421) 

— astrometry 


1. INTRODUCTION 

Blazars are extremely bright active galactic nuclei 
(AGNs) with variable jets. Blazar jets are thought to be 
powered by the accretion of surrounding material onto 
super-massive black holes (“central engines”). The radio 
jet bases (“radio cores”) located in the upstream ends 
of AGN jets in radio images have been basically consid¬ 
ered to be opaque photospheres (surfaces of an optical 
depth of ^ 1) against synchrotron self-absorption (SSA) 
or free-free absorption. Frequency-dependent offsets of 
the radio core position detected by very long baseline 
interferometry (VLBI) observations support this charac¬ 
teristic structure (e.g., Marcaide & Shapiro 1983; Hada 
et al. 2011). On the other hand, the stationary stand¬ 
ing shock model, which is also another idea explaining 
blazar radio cores has been proposed based on recent 
VLBI observations (e.g., Marscher et al. 2008). Although 
detailed radio observations can provide the morpholog¬ 
ical structure of an inner jet, their stationarity is less 
well understood in addition to the location of the radio 
core still being a matter of debate (e.g., Marscher et al. 
2008; Hada et al. 2011). Long-period VLBI astrometry 
at low frequency has not only derived the stationarity of 
radio cores with an astrometric accuracy (e.g., Bartel et 
al. 1986; Guirado et al. 1995; Fey et al. 1997) but also 
caught the sign of their wandering behavior (Ros et al. 
1999; Marti-Vidal et al. 2011). In particular, the sta¬ 
tionarity of blazars related to high-energy activity has 
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not yet been verified by observations within a limited 
period after the flaring phenomena. 

The high-synchrotron-peaked BL Lac object Markar¬ 
ian 421, at redshift z = 0.031, is known as the first bright 
7 -ray blazar in the TeV energy range to be reported by 
the Whipple ground-based Gherenkov telescope (Punch 
et al. 1992). This blazar harbors a ~ 3.6 x 10® solar mass 
black hole (Wagner 2008). Niinuma et al. (2012) have re¬ 
ported on superluminal apparent inward motion of the 
jet component relative to the radio core, which has steep 
spectra at more than ~10 GHz (Lico et al. 2012; Blasi 
et al. 2013), using the Japanese VLBI network (JVN) 
soon after the largest X-ray flare occurred in 2010. It is 
clear that inward motions along outflowing jets are not 
realistic, but these apparent motions can be explained if 
the position of the radio core changes (Kellermann et al. 
2004). 

To clarify the origin of such inward motions, we per¬ 
formed 12 epochs of astrometric observations of the large 
X-ray flare of Markarian 421 in 2011 September (Hiroi 
et al. 2011) using the VLBI Exploration of Radio As¬ 
trometry (VERA) at 22 GHz together with a nearby 
reference source JllOl-l-3904 (see Fig. 1). Through¬ 
out this paper, we adopt cosmological parameters of 
Ho = 70.2 kms“^ Mpc“^, Da = 0.725 (Komatsu et al. 
2011). Based on these parameters, 1 mas at z = 0.031 
corresponds to a linear scale of 0.62 pc. 

2. ASTROMETRIG OBSERVATION AND DATA ANALYSIS 
2.1. VERA observations 

We conducted 12 epochs of phase-referencing observa¬ 
tions at 22GHz between 2011 September 16 and 2012 
April 28 with typical interval of 20days after the afore¬ 
mentioned large X-ray flare. Observations were con¬ 
ducted using VERA (baselines ranging between 1000 and 
2300 km). VERA was designed for dedicated astrometric 
observations by installing a dual beam system (beams A, 
B) and is operated by the National Astronomical Obser¬ 
vatory of Japan. 

To measure the absolute position of Markarian 421, we 
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Fig. 1.— Configuration of Markarian 421 and JllOl+3904 jets 
on the celestial sphere. The tracking centers for both Markar- 
ian 421 and J1104+3904 are R.A. = ll‘'04™27=.313943, decl. 
= +38° 12'31".79906, and R.A. = llhoi"°30+069578, decl. = 
+39°04'32”.63332, respectively. The separation angle of both 
sources is 1.04°. Red arrows represent jet orientations of both 
sources. 

selected JllOl+3904, which is listed in the VLBA cali¬ 
brator survey list^, as a position reference. In Fig. 1, we 
represent the configuration of both sources on the celes¬ 
tial sphere. To clearly display the source morphologies 
and jet orientations of Markarian 421 and JllOl-l-3904, 
we show 8 GHz VLBI images of both sources observed 
by the Very Long Baseline Array (VLBA)®. The position 
angle (P-A.) measured from north to west of the jet of 
JllOl-l-3904 is 153°and that of Markarian 421 is 35°on 
average, respectively. These jet orientations were derived 
from VERA images and are indicated with red arrows in 
Fig. 1. Therefore, the two jet axes are nearly perpendic¬ 
ular (118°). In addition, the separation angle between 
two sources on the celestial sphere is only 1°.04. The or¬ 
thogonality and the closeness of two sources are crucial 
to avoid any degeneracy in the structural effects in the 
jet direction and to reduce the astrometric errors caused 
by the tropospheric delay residuals. Therefore, the selec¬ 
tion of this pair of jets allows us to measure the absolute 
position of Markarian 421 with respect to JllOl-l-3904 
with high astrometric accuracy. 

Our VERA observations were performed with a record¬ 
ing rate of 1024Mbps (2-bit quantization x 16 MHz x 16 
intermediate frequency (IF) channels ). Eight IF chan¬ 
nels were assigned to Markarian 421 (beam A), and the 
others were assigned to JllOl-l-3904 (beam B). The in¬ 
strumental phase difference between the two beams was 
monitored in real time during the observations by in¬ 
jecting artificial noise sources into both beams at each 
station (Honma et al. 2008a). 

2.2. Data analysis 

The visibility calibration was performed using the 
NRAO Astronomical Image Processing System (AIPS) 
software package. To increase astrometric precision, we 
recalculated the delays of all sources using a precise 
geodetic model and calibrated the instrumental phase 
difference between the dual beam systems. This precise 

^ http://www.vlba.nrao.edu/astro/calib/vlbaCalib.txt 
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Fig. 2.— Left and right panels show, respectively, the self- 
calibrated and phase-referenced VERA images of J1101-h3904 on 
2012 April 6. These images are convolved with a circular Gaus¬ 
sian function that has a full-width at half-maximum of 0.8 mas, the 
typical minor-axis size of the synthesized beam in our observations, 
shown in the bottom-left corner. 

recalculation of the delay also included the most recent 
Earth rotation parameters provided by the International 
Earth Rotation and Reference System Service (lERS) 
and the tropospheric delays measured using global posi¬ 
tioning system receivers at each VERA station (Honma 
et al. 2008b). Ionospheric delays were considered based 
on a global ionosphere map, which was produced every 
2hr by the University of Bern. 

We performed fringe-fitting on both Markarian 421 and 
JllOl-l-3904 independently using the AIPS task FRING 
and removed the residual delays, rates, and phases as¬ 
suming a point-source model. Then VLBI images of 
Markarian 421 and JllOl-l-3904 (Fig. 2 left) were pro¬ 
duced via a standard GLEAN and self-calibration proce¬ 
dure using the Caltech Difmap package (Shepherd 1997). 

To obtain the core position of Markarian 421, we 
performed a phase-referencing analysis in the following 
manner. First, we chose Markarian 421 as the phase 
calibrator because this source is much brighter than 
JllOl-l-3904. After the a priori corrections of the visi¬ 
bility amplitude, the solutions for residual delays, rates, 
and phases obtained on Markarian 421 (under the point- 
source assumption) were transferred to the data set of 
JllOl-l-3904. Then, to account for the source structure 
phase of Markarian 421, we ran the AIPS task CALIB on 
Markarian 421 with the source model derived from the 
self-calibrated image, and the solution was also trans¬ 
ferred to JllOl-l-3904 for accurate determination of the 
absolute position of the radio core. Based on these pro¬ 
cedures, we finally obtained phase-referenced images of 
JllOl-l-3904 (Fig. 2 right), which displayed the relative 
position of JllOl-l-3904 with respect to Markarian 421. 
In the astrometry result, we present the relative position 
of Markarian 421 with respect to JllOl-l-3904. 

To determine the jet orientation of both sources uti¬ 
lized through this paper, we fitted two-dimensional cir¬ 
cular Gaussian models to these self-calibrated visibilities. 
The jet orientation of Markarian 421 derived from the 
VERA observations is consistent with that reported by 
previous VLBA observation (Lico et al. 2012). The P.A. 
of the innermost jet in JllOl-l-3904 (~0.8 mas from the 
core) derived from VLBA 8 GHz data (see Fig. 1) is 
approximately 150°, which is also quite consistent with 
VERA observations. 

The astrometric accuracy in each observing epoch rep¬ 
resented in Fig. 3-(b) and 4 is 0.08-0.12 mas. To esti- 
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Fig. 3.— Measured radio-core position of Markarian 421 with re¬ 
spect to the position of the reference source J1101-h3904. (a): The 
horizontal and vertical axes and the color gradient represent the 
positional offsets of right ascension and declination from the obser¬ 
vation coordinates with one-sigma positional uncertainties and the 
observation date, respectively. The black arrow indicates the jet 
direction of Markarian 421. (b): The horizontal and vertical axes 
show the time and the positional offsets along the perpendicular 
direction to the jet of J1101-1-3904. The black arrow indicates the 
date on which the large X-ray flare occurred. 

mate these values, we calculated (1) the statistical astro- 
metric error derived from the signal-to-noise ratio of the 
phase-referenced image of J1101-1-3904, (2) the dispersive 
ionospheric residuals, (3) the non-dispersive tropospheric 
residuals, and (4) the identification of the core position 
by following the way described in Hada et al. (2011). 

3. DETECTION OF THE POSITION WANDERING OF THE 
RADIO JET BASE 

We discovered a core-position wandering in Markar¬ 
ian 421 based on bonafide astrometric observations (Fig. 
3). This is the first direct detection of such a clear po¬ 
sition change of the radio core within a year at higher 
frequency (> 20 GHz). We show the position of the 
Markarian 421 radio core in the celestial coordinates in 
Fig. 3-(a). The core clearly shows its position change 
in the northwest direction, which is a similar direction 
to the Markarian 421 jet axis. In Fig. 3-(b), we present 
the measured core positions of Markarian 421 projected 


to the P.A. of 63°, which is a direction perpendicular to 
the JllOl-l-3904 jet axis, because the structural effect of 
J1101-1-3904 is negligible in the motion along this direc¬ 
tion. The position of the radio core began to change sig¬ 
nificantly ^ 20 days after the peak of the X-ray flare. The 
radio core, which is located at —0.16 ±0.11 mas on 2011 
September 16 (modified Julian date (MJD) 55820) in Fig. 
3-(b), moved in the positive direction to 0.35 ± 0.11 mas 
by 2011 October 25 (MJD 55859). The core position 
changed by 0.51 mas within 39days, which corresponds 
to an apparent velocity of 9 c. After the core returned 
to near the initial position on MJD 55890, a core-position 
change of nearly the same magnitude, 0.49 mas, was de¬ 
tected again within the 41 days between MJD 55890 and 
MJD 55931. The radio core then gradually returned once 
more to its initial position by approximately MJD 56000. 
We stress that the amount of core-position wandering 
is significantly larger than our astrometric accuracy of 
^ 100 fj,as. 

For quantitative evaluation of the significance of the 
radio core wandering during all the observing epochs, we 
examine the chi-square test for the position change rep¬ 
resented in Fig. 3-(b). In this test, we define two cases of 
null hypothesis as follows: the position of the radio core 
shows the stationarity for 1) the weighted mean position 
of the radio core during all the observing epochs and 2) 
the weighted mean position of the radio core during last 
three epochs in which it seems to be stationary. In both 
cases, the stationarity of the radio core is rejected at a 
significance level of < 0.1% as shown in Table 1. The 
radio core, hence, exhibits its position change during ob¬ 
serving epochs at a confidence level of > 99.9%. 

In Fig. 4, we present the absolute core positions along 
the jet axis of Markarian 421, and radio images are over¬ 
laid on each corresponding epoch. In this figure, the 
presence of a single jet component (hereafter, we call it 
‘Cl’) located at approximately 1.3 mas northwest from 
the core can be confirmed and is also well known as its 
stationarity (Lico et al. 2012). Cl is much fainter than 
the radio core as represented in Fig. 4. In three out of 
twelve epochs (namely 2011 October 25, 2011 December 
16, and 2012 January 05), the radio core significantly 
moves downstream (in the northwest direction) of the 
jet, and the separations between the bright core and Cl 
become too close to resolve these two components with 
VERA’s synthesized beam. Hence, Cl was not seen in 
these three epochs. Thus, the 12 epoch sequence of the 
radio images suggests the core position wandering phe¬ 
nomenon discovered by our astrometric observations. 

Moreover, if a displacement of the brightness centroid 
between the stationary core and a newborn jet causes 
the core wandering phenomenon, an elongated or double- 
peaked radio core image around the reference stationary 
position (the white line), which should be expected in the 
epochs where the core clearly moves downstream of the 
jet in Fig. 4 (e.g., radio images with a beam size of 0.1 
mas in Marscher et al. 2010). However, such signatures 
were not observed in the radio-core images. Therefore, 
the wandering behavior of the radio core discovered in 
the present study clearly differs from the observational 
feature of the radio core in, e.g., BL Lacertae, in which 
the radio core is indicated to be stationary (Marscher et 
al. 2008). 
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Fig. 4. — Radio-core positions projected on the jet axis and radio images of Markarian 421. These images are overlaid on their corre¬ 
sponding absolute core positions and were convolved with a circular Gaussian function that has a full width at half maximum of 0.8 mas 
(shown in the bottom-right corner). The contours begin at 6 mjy beam“^ and increase in steps of 2"/^. The vertical axis represents the jet 
axis. The white arrow in 2011 September indicates the date of the large X-ray flare. The white solid and the blue dashed lines indicate the 
positions of the radio core and the jet component (Cl) as measured in 2011 September 16, respectively. The yellow dashed line represents 
the cubic spline interpolation of each core position. 


TABLE 1 

Chi-square test for stationarity of the radio core 

POSITION. 


Parameters 

Case (1) 

Case (2) 


0.002 

-0.106 


40.99/11 

56.00/11 

-^stationary (X^ 5 

< 0.1 

< 0.1 


Note. 

t/r is the weighted average position of the radio core projected to 
the P.A. of 63°in mas. 
tu is the number of degrees of freedom. 

Probability (%) that the radio core is stationary for the position 

4. DISCUSSION 

The de-projected distance that corresponds to the 
core-position changes of Markarian 421 is approximately 
3.5 ^ 8.7 pc (equivalent to (1.0 — 2.6) x 10^ Schwarzschild 
radii (Rs)) for a jet-viewing angle of 2 — 5°(Lico et al. 
2012). One of the leading models for explaining the 
blazar phenomenon is the so-called internal-shock model, 
in which discrete ejecta with higher speeds catch up with 
the preceding slower ejecta, and the collision leads to a 
shock wave in the colliding ejecta (Spada et al. 2001; 
Guetta et al. 2004). Based on this model, core wander¬ 
ing can be naturally explained because the model pre¬ 
dicts some scattering of the shock locations caused by 
the randomness of the ejecta speed. Therefore, we as¬ 


sume that the distances of the shock locations from the 
black hole immediately after the flare, D ^ 10® i?s, are 
proportional to the product of the square of the bulk 
Lorentz factor (Tjet) of the ejecta and the separation 
length (A) of the colliding ejecta. Assuming that a typi¬ 
cal separation length is a scale of about 10 times the in¬ 
nermost stable circular orbit for non-rotating black holes 
(fisco = 3 i?s), in light of numerical simulations of rel¬ 
ativistic outflows around black holes (McKinney et al. 
2013), very fast ejecta with a Lorentz factor of approx¬ 
imately Ljet 60 (H/IO® i?s)^^^(A/30 are re¬ 

quired to explain the observed magnitude of the core 
wandering. The Tjet is comparable to previously reported 
maximum values (T of 40 ^ 50; Lister et al. 2009, 2013) 
based on statistical studies by a long-term monitor using 
VLBI. Moreover, the radio core returned to the reference 
stationary position after 2012 March 12 (the 10th epoch). 
This can be explained if the emission from the reference 
position is identical to the radio core that is in a quies¬ 
cent state. The radio-core emission in a quiescent state 
is expected to be persistent. As one of the possibilities 
for reproducing the observed core wandering, this persis¬ 
tent radio core was obscured by the flare-associated radio 
core, which therefore must be optically thick against SSA 
and short lived (a lifetime of approximately 2 weeks), for 
several epochs. Such a flare-associated radio core must 
be located at least 10® Rs downstream from the per¬ 
sistent radio core. In addition, Koyama et al. (2015) 
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recently conducted multi-epoch astrometric observations 
of Markarian 501, which is also one of the most nearby 
TeV blazars by using VERA. However, they detected no 
significant positional changes of the radio core during its 
quiescent state. This fact, therefore, seems to support 
the radio core wandering seen after the large X-ray flare. 

Another possible scenario, which may lead to the radio¬ 
core wandering is the presence of changes in the elec¬ 
tron density and magnetic field strength in the jet flow. 
Such changes might be triggered when magnetohydrody¬ 
namic waves propagate in the jet (Rees 1971; Cohen et al. 
2014) with an apparent wave-pattern speed of approxi¬ 
mately 9 c. The propagation of electromagnetic waves in 
a charge-starved jet proposed by Kirk & Mochol (2011) 
might also be generated at the launching point of the jet, 
although an assumption of charge starvation may not be 
warranted because of the high number density of low- 
energy electrons (Kino & Takahara 2008). In addition, 
Gomez et al. (1997) have predicted that the interaction 
between the steady standing shock and a propagating 
perturbation that contains shocked regions and is trailed 
by a rarefaction along the jet causes a positional shift of 
the brightness centroid with a superluminal speed com¬ 
parable to the one reported in this Letter. 

Regardless of the core-wandering mechanism, it is ev¬ 
ident that our new observational finding offers crucial 
hints to resolving a key jet origin question (i.e., particle 
versus electromagnetic power) in the innermost region of 
the TeV blazar Markarian 421. In addition, recent radio 
and optical observations have suggested that monitoring 
the behavior of polarized emission on a short timescale 
is also a promising way to understand the origin of the 
jet activity in the innermost region of blazars (Lico et al. 
2014; Itoh et al. 2015). 

Thus far, as one of the few examples, multi-epoch as¬ 
trometric observations have been performed for the FR- 
1-type radio galaxy M87 (Acciari et al. 2009; Hada et al. 
2014). Acciari et al. (2009) noted that the position of 
the M87 radio core at 43GHz moved by no more than 
6 Rs during the flare that occurred in 2008, while Hada 
et al. (2014) similarly showed a stable (less than 40 Rs) 
nature for the 43GHz core position during the flare in 
2012. By contrast, we report here the detection of a large 
change in the position of a radio core. The de-projected 
distance corresponding to the core-position changes of 
Markarian 421, (1.0 — 2.6) x 10® i?s, is approximately 
ten thousand times larger than the upper limit in the 
case of M87. According to the AGN unified model (Urry 
& Padovani 1995), both the blazar Markarian 421 and 
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the radio galaxy M87 are classified as FR-l-type galax¬ 
ies, and their jet power is approximately of Pjet ^ 10®^ 
W (Rieger & Aharonian 2012; Potter & Cotter 2013). 
Therefore, the difference is most likely related to the jet¬ 
viewing angle and the optical depth structure rather than 
the jet power. The results reported by Piner et al. (2010) 
and Hada et al. (2013) imply that the jet has a transverse 
structure composed of an inner fast spine flow and a sur¬ 
rounding slow layer (Ghisellini et al. 2005). On the basis 
of this model, observers see the radio core in the blazar, 
which is dominated by the fast spine flow because the 
jet is closely aligned with our line of sight. In contrast, 
the radio galaxy has a jet off-axis whose emission comes 
from the slow layer and not from inner spine flow. A 
structured jet is a possible model for explaining the ob¬ 
servational discrepancy in the behavior of the radio cores 
between Markarian 421 and M87. 

5. SUMMARY 

We conclude that the radio core-wandering phe¬ 
nomenon implies that the blazar emission zone 
in Markarian 421 is occasionally located at least 
(1.0 — 2.6) X 10® Rs downstream from the black hole just 
after the large X-ray flare. To wrap up, the location of 
the radio core in Markarian 421 is far from the black hole, 
as also suggested in other blazars but without perform¬ 
ing astrometric observations (e.g., Marscher et al. 2008, 
2010; Abdo et al. 2010; Orient! et al. 2013), while the 
radio core is not stationary after the large X-ray flare. 
In addition, our result on the location of the radio core 
is consistent with the one in the case of GTA 102, which 
was classified as a flat spectrum radio quasar measured 
by the image alignment method (Fromm et al. 2015, and 
references therein). This radio core wandering was dis¬ 
covered just after the large X-ray flare phenomenon us¬ 
ing dense VLBI astrometric observations. Thus there is 
a possibility that the wandering behavior may be a new 
type of manifestation associated with the phenomena of 
large X-ray flares. 
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